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Abstract. A hypothetical time-variation of the gravitational constant G would make neutron 
stars expand or contract, so the matter in their interiors would depart from beta equilibrium. 
This induces non-equilibrium weak reactions, which release energy that is invested partly in 
neutrino emission and partly in internal heating. Eventually, the star arrives at a stationary 
state in which the temperature remains nearly constant, as the forcing through the change 
of G is balanced by the ongoing reactions. Using the surface temperature of the nearest 
millisecond pulsar (PSR J0437-4715) inferred from ultraviolet observations and results 
from theoretical modelling of the thermal evolution, we estimate two upper limits for this 
variation: (1) \G/G\ < 2 x 10~ 10 yr -1 , if the fast, "direct Urea" reactions are allowed, and (2) 
\G/G\ < 4 x 10~ 12 yr -1 , considering only the slower, "modified Urea" reactions. The latter 
is among the most restrictive upper limits obtained by other methods. 
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1. Introduction 

A number of theorists have proposed that the 
so-called "fundamental constants" of Nature 
may vary with cosmological time, ever since 
Dirac suggested th at the gra vitational force 
may be weakening dDiradll937h . Several ex- 
periments aimed at constraining the time vari- 
ation of G have been conducted, and their 
results can be grouped by the time scales 
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they probe (see iReisenegger et al.|[2007l for a 
list with references): (1) Human lifetime (~ 
10 yr) experiments rely on real-time moni- 
toring of distances within the Solar System, 
white-dwarf oscillation periods, and pulse ar- 
rival times of isolated and binary pulsars, (2) 
long time (~ 10 9 ~ 10 yr) measurements use stel- 
lar astrophysics and paleontology, and (3) cos- 
mological (> 10 10 yr) experiments use the Big 
Bang nucleosynthesis and Cosmic Microwave 
Background anisotropics to compare the value 
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of G in the early Universe to the present value. 
Generically, they set constraints on \G/G\ down 
to ~ 1(T 12 yr _1 , although the comparison be- 
tween different timescales depends on the as- 
sumed form of the function G(t). 

Here, we review our previously introduced 
metho d for setting constraints on G dJofre et al.l 
2006), to our knowledge the only one so far 
to probe timescales ~ 10 7 ~ 9 yr, intermediate 
between the human and long timescales men- 
tioned above. It relies on the change in the in- 
ternal structure of a neutron star induced by a 
change in G, which, together with the slow re- 
sponse timescale of weak interactions, result in 
internal heating and an increase in the observed 
surface temperature. 

Neutron stars have mean densities exceed- 
ing nuclear saturation density, p nuc ~ 3 x 
10 14 g cm 4 . In their outer layers, they are com- 
posed of heavy atomic nuclei and free elec- 
trons, giving way to free neutrons, free protons, 
muons, and potentially more exotic particles as 
the density increases inward. A short time af- 
ter their formation, their internal temperatures 
drop orders of magnitude below the Fermi en- 
ergies of free particles (~ 10 - 100 MeV), and 
so their structure is well approximated by zero- 
temperature models. 

Nonetheless, weak interactions still play an 
important role, as characteristic temperatures 
10 6 ~ 8 K yield non-negligible neutrino emis- 
sion. Given that the equation of state of mat- 
ter above nuclear density still remains poorly 
known, researchers construct models to cal- 
culate the evolution of the thermal content of 
neutron stars and compare their predictions 
with X-ray observations, in order to set con- 
straints on the models for dense matter (e.g., 
lYakovlev & Pethickl |2004|) . In these thermal 
evolution models, the structure of the star is 
calculated assuming an equation of state, and 
the evolution of the internal temperature is ob- 
tained by considering losses due to different 
neutrino emission processes from the stellar in- 
terior, as well as thermal electromagnetic radi- 
ation from the surface. 

We have previously explored rotochemi- 
cal heating, the effect of the progressive loss 
of rotational support of millisecond pulsars 
has on their internal structure, which leads 



to weak interactio n processes (beta d ecays) 
and net heating (Reisenegger] 1 19951 Il997t 
iFernandez & Re isenegger 2005). A time varia- 
tion of the gravitational constant has an analo- 
gous effect: it compresses the star and therefore 
also causes heat ing (gravitochemical heating, 
IJofre et al.ll2006h . In what follows, we discuss 
how this process constrains G. 



2. Method 

Neutron star matter is composed of degenerate 
fermions of various kinds: neutrons (n), pro- 
tons (p), leptons (/), and possibly other, more 
exotic particles. Neutrons are stabilized by the 
presence of other, stable fermions that block 
(through the Pauli exclusion principle) most 
of the final states of the beta decay reaction 
n — > p + I + v, making it much slower than 
in vacuum. The large chemical potentials pi 
Fermi energies) for all particle species i also 
make inverse beta decays, p + /—»« + v, pos- 
sible. The neutrinos (v) and antineutrinos (v) 
leave the star without further interactions, con- 
tributing to its cooling. The two reactions men- 
tioned tend to drive the matter into a chemical 
equilibrium state, defined by r] np i = p n - p p - 
pi = 0. Depending on the relative abundance of 
protons, the reactions mentioned above (called 
direct Urea processes) may not be able to 
conserve momentum, in which case an addi- 
tional nucleon must participate in the reactions 
(yielding the modifie d Urea processes, e. g., 
lYakovlev et al.l 1200 ll) . with strongly reduced 
rates. 

If G changes, so does the hydrostatic equi- 
librium structure of cold neutron stars, and 
all matter elements in their interior are com- 
pressed or expanded, and in this way driven 
away from chemical equilibrium (r\„ p i + 0). 
Free energy is stored, which is released by an 
excess rate of one reaction over the other. This 
energy is partly lost to neutrinos and antineu- 
trinos (undetectable at present), the remainder 
heats the stellar interior. The heat is eventually 
lost as thermal (typically ultraviolet) photons 
emitted from the stellar surface. 

The thermal evolution of the star is calcu- 
lated by solving a system of ordinary differen- 
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tial equations of the form 
t = -J^ [l h (T, r] np i) - L V (T, r) np i) 

~ Ly{T)\ 



C(T) 

f/npi = A{T,t] np i) + BG 



(1) 

(2) 



where T is the internal temperature, C the heat 
capacity of the star, L#, L v , and L y the heating, 
neutrino cooling, and photon cooling luminosi- 
ties, respectively, A represents the change in 
the particle abundances due to reactions, and 
B a constant coefficient tha t depends only on 
the structure of the star. See lJofre et a (120061) 
and lFernandez & Rei senegger (2005|)f6r more 
details. 



3. Results 

The typical result of integrating equations ([TJ 
and (f2]i is as follows (details for the 
analogou s, rotochemical heating case are 
given in [Fernandez & Reisenegger! 120051 and 
iReisenegger et al.l 2007 ). The star first cools 
down (within < 10 7 yr) from its high birth 
temperature, while the chemical potential im- 
balances rj np i slowly increase due to the grav- 
itational forcing term in equation (fJJ. After 
~ 10 7 ~ 9 yr, depending on \G\, the imbalances 
are so high that they increase the reaction rates 
to the point where they keep up with the forc- 
ing, stabilizing the chemical imbalance. At this 
point, the energy released by reactions and that 
emitted from the surface also compensate in 
equation ([TJ, keeping the temperature at a con- 
stant value ~ 10 5 K, whose precise value can 
be predicted if \G\ and the neutron star model 
(mass and equation of state) are given. 

In order to constrain G, we therefore need 
sensitive ultraviolet observations of known 
neutron stars older than 10 7 yr. Such an obser- 
vation only exists for the nearest millisecond 
pulsar, PSR J0437-4715, for which an HST 
observation detected what appears to be the 
ultraviolet Rayleigh- Jeans tail of a blackbody 
at ~ 10 5 K dKargaltsev et al.ll2004 . Figures Q] 
and [2] compare theoretical predictions for var- 
ious equations of state and a range of neutron 
star masses to the error contours obtained from 
this observation. The upper limit on \G/G\ = 
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Fig. 1. Comparison of the gravitochemi- 
cal heating predictions to observations of 
PSR J0437-4715. The solid lines are the pre- 
dicted stationary surface temperatures as func- 
tions of stellar radius, f or different equa- 
tions of state ("A 18 from lAkmal et al.l 1 19981 



and BPAL from iPrakash et alj 1 19881) . con- 
strained to the obse rvationally allowed mas s 
range for this pulsar Ivan Straten et al.l d200 lh . 
Dashed lines correspond to the 68% and 
90% confidence contours o f the blackbody 
fit of fkargaltse v et alJ (120041) for the ultravi- 
olet emission from this object. The value of 
\G/G\ = 2x 10~ 10 yr~' is chosen so that all sta- 
tionary temperature curves lie above the obser- 
vational constraints. (BPAL32 and BPAL33 al- 
low direct Urea reactions in the observed mass 
range of J0437.) 
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Fig. 2. Same as Figure [TJ but now the value of 
\G/G\ - 4 x 10~ 12 yr -1 is chosen so that only 
the stationary temperature curves with modi- 
fied Urea reactions are above the observational 
constraints. 
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2 x 1CT 10 yr -1 (Fig. [0 is obtained by requiring 
that all stationary temperature curves lie above 
the observational constraints. Restricting the 
models to only modified Urea reactions yields 
a more restrictive constraint, \G/G\ = 4x 10" 12 
yr" 1 (Fig.0. 
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Yakovlev, D. G., et al. 2001, ApJ, 354, 1 
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4. Discussion 

Gravito chemical heating sets constraints on 
the time variation of the gravitational con- 
stant, using the fact that a non-zero change 
would generate internal heating in neutron 
stars, which for nearby cases can be detected 
with existing telescopes. These constraints are 
the only ones on timescales 10 7 ~ 9 yr. If it 
could be assured that the observed neutron 
stars cannot have direct Urea reactions, these 
constraints would be of the same order as the 
best ones available on other time scales. This is 
not the case at the moment, but could become 
so as neutron star interiors become better un- 
derstood from other studies. In particular, ob- 
served temperatures of other old neutron stars 
(such as millisecond pulsars) will provide use- 
ful information. 
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